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A general formalism is set up to analyse the response of an arbitrary solid elastic body to an 
arbitrary metric Gravitational Wave perturbation, which fully displays the details of the interaction 
antenna-wave. The formalism is applied to the spherical detector, whose sensitivity parameters are 
thereby scrutinised. A multimode transfer function is defined to study the amplitude sensitivity, and 
absorption cross sections are calculated for a general metric theory of GW physics. Their scaling 
properties are shown to be independent of the underlying theory, with interesting consequences for 
future detector design. The GW incidence direction deconvolution problem is also discussed, always 
■ within the context of a general metric theory of the gravitational field. 
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I. INTRODUCTION 



The idea of building ultracryogenic spherical Gravitational Wave (GW) antennae seems to be progressively winning 
adepts, even despite the technological difficulties of various kinds posed by a project like that, which every expert 
■ acknowledges. Confidence in its feasibility stems from many years of experience: groups at Stanford, Louisiana State 
University, Roma and Legnaro (Italy) and Perth (Western Australia) have constructed and operated, at different 
levels, cryogenic cylindrical bars of the Weber type In particular, a long term strain sensitivity h — 6 x 10~ 19 
for millisecond bursts has been reported from the bar EXPLORER Q. The new generation ultracryogenic cylinder 
NAUTILUS, of the Frascati group jj), is beginning operation as these lines are written ||, with an expected sensitivity 
nearly an order of magnitude better than the above. 

Spherical antennae are considered by many to be the natural next step in the development of resonant GW detectors 
|)| |l0|. The reasons for this new trend essentially derive from the improved sensitivity of a sphere — which can be 
nearly an order of magnitude better than a clyinder having the same resonance frequency, see below and |h| — , and 
?H . from its multimode capabilities, first recognised by Forward (j| and further elaborated in |@,||. 

Although some of the most relevant aspects of detector sensitivity have already received attention in the literature, it 
J> , seems to me that a sufficiently general and flexible analysis of the interaction detector-GW has not been satisfactorily 
developed to date. This theoretical shortage has a number of practical negative consequences, too. Traditional analysis, 
to mention but an example, is almost invariably restricted to General Relativity or scalar-tensor theories of gravity; 
while it may be argued that this is already very general, any such argument is, as a matter of fact, understating the 
potentialities actually offered by a spherical GW antenna to help decide for or against any one specific theory of the 
gravitational field on the basis of experimental observation. 

I thus propose to develop in this paper a full fledged mathematical formalism which will enable analysis of the 
antenna's response to a completely general GW, i.e., making no a priori assumptions about which is the correct 
theory underlying GW physics (other than, indeed, that it is a metric theory), and also making no assumptions about 
detector shape, structure or boundary conditions. Considering things in such generality is not only "theoretically 
nice" — it also brings about new results and a better understanding of older ones. For example, it will be proved 
that the sphere is the most efficient GW elastic detector shape, and that higher mode absorption cross sections scale 
independently of GW physics. I will also discuss the direction of incidence deconvolution problem in the context of a 
general metric theory of gravity. 

The paper is organised as follows: section 2 is devoted to the development of the general mathematical framework, 
leading to a formula in which an elastic solid's response is related to the action of an arbitrary metric GW impinging 
on it. In section 3 the general equations are applied to the homogeneous spherical body, and a discussion of the 
deconvolution problem is presented as well. Section 4 contains the description of the sphere's sensitivity parameters, 
specifically leading to the concept of multimode, or vector, transfer function, and to an analysis of the absorption 
cross section presented by this detector to a passing by GW. Conclusions and prospects are summarised in section 5, 
and two appendices are added which include mathematical derivations. 
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II. GENERAL MATHEMATICAL FRAMEWORK 



In the mathematical model, I shall be assuming that the antenna is a solid elastic body which responds to GW 
perturbations according to the equations of classical non-relativistic linear Elasticity Theory [jfl). This is fully justified 
since, as stressed above, GW induced displacements will be very small indeed, and the speed of such displacements 
much smaller than that of light for any forseeable frequencies. Although our primary interest is a spherical antenna, 
the considerations which follow in the remainder of this section have general validity for arbitrarily shaped isotropic 
elastic solids. 

Let u(x, t) be the displacement vector of the infinitesimal mass element sitting at point x relative to the solid's 
centre of mass in its unperturbed state, whose density distribution in that state is p(x). Let A and p, be the material's 
elastic Lame coefficients. If a volume force density f(x, t) acts on such solid, the displacement field u(x, t) is the 
solution to the system of partial differential equations jllj 

p|^MV 2 u-(A + M )V(V-uHf(x,i) ( 2J ) 

with the appropriate initial and boundary conditions. A summary of notation and general results regarding the solu- 
tion to that system is briefly outlined in the ensuing subsection, as they are necessary for the subsequent developments 
in this paper, and also in future work. 



A. Separable driving force 

For reasons which will become clear later on, we shall only be interested in driving forces of the separable type 

f(x,t) = f(x) ff (t) (2.2) 



or, indeed, linear combinations thereof. The solution to (2.1) does not require us to specify the precise boundary 
conditions on u(x, t) at this stage, but we need to set the initial conditions. We adopt the following: 

u(x,0) = u(x,0) = (2.3) 

where '= d/dt, implying that the antenna is at complete rest before observation begins at t=0. The structure of 
the force field (2.2) is such that the displacements u(x, t) can be expressed by means of a Green function integral of 
the form 



/>oo 

u(x,i) = / S(x;t-t')g(t')dt' 
Jo 



(2.4) 



The deductive procedure whereby S(x; t — t') is calculated can be found in many standard textbooks — see e.g. |l2| . 
The result is 

r if t < 

S(x;t)=i (2.5) 

[ Eat ^ Uat(x) Sin^Art if t > 

where 

fN=lf I u^(x).f(x)d 3 x (2.6) 

m J Solid 

and utv(x) are the normalised eigen-solutions to 

^V 2 uat + (A + a*)V(V-ujv) = -u N pu N (2.7) 

with suitable boundary conditions. Here N represents an index, or set of indices, labelling the eigenmode of 
frequency lun- The normalisation condition is (arbitrarily) chosen so that 

/ u^,(x) • u N (x)p(x)d 3 x = MS N ' N (2.8) 

./Solid 
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wh ere M is the total mass of the solid, and the asterisk denotes complex conjugation. Replacing now ( |2.5| ) into 
( |2.4|) we can write the solution to our problem as a series expansion: 

u(M)=£— Mx) flw (t) (2.9) 



where 



g N (t) = [ g{t') sin uj n (t-t') dt' (2.10) 



Equation (|]^) is the formal solution to our proble m; it has the standard form of an orthogonal expansion and 



is valid for any solid driven by a separable force like (2.2) and any boundary conditions. It is therefore completely 
general, given that type of force. 

Before we go on, it is perhaps interesting to quote a simple but useful example. It is the case of a solid hit by a 
hammer blow, i.e., receiving a sudden stroke at a point on its surface. Exam of the response of a GW antenna to 
such perturbation is being used for correct tuning and monitoring of the device If the driving force density is 
represented by the simple model 

f( hb >(x,i)=fo^(x-x )<5(t) (2.11) 
where Xo is the surface point hit, and fo is a constant vector, then the system's response is immediately seen to be 

u (hb) (x,i) =V ^-ujv(x) sin cj N t (2.12) 

N " 

with fff = M^ 1 fo-u^(xo). A hammer blow thus excites all the solid's normal modes, except those perpendicular 
to fo , with amplitudes which are inversely proportional to the mode 's frequency. This is seen to be a rather general 
result in the theory of sound waves in isotropic elastic solids. 

B. The GW tidal forces 



An incoming GW manifests itself as a tidal force density; in the long wavelength linear approximation 14 it only 
depends on the "electric" components of the Riemann tensor: 

fi(x,t)=pc 2 R oi oj(t)x j (2.13) 



where c is the speed of light, and sum over the repeated index j is understood. In ( [2.15 ) tidal forces are referred 



to the antenna's centre of mass, and thus x is a vector originating there. Note that I have omitted any dependence 
of Roioj on spatial coordinates, since it only needs to be evaluated at the solid's centre. The Riemann tensor is only 
required to first order at this stage Jl5| : 

RoiOj — 2 {hijfiO — hoifij — hojfii + hoo,ij) (2-14) 
where h^ v are t he perturbations to flat geometry^, alway s at the centre of mass of the detector 



The form ( 2.13 ) is seen to be a sum of three terms like (2.2) — but this three term "straightforward" splitting is 
not the most convenient, due to lack of invariance and symmetry. A better choice is now outlined. 
An arbitrary symmetric tensor Sij admits the following decomposition: 

2 

S ij (t)=SV\t)E$P + S (m) (t)4f ( 2 - 15 ) 

m=— 2 



1 Throughout this paper, greek indices (/i, v, . . .) will run through space-time values 0,1,2,3; latin indices . . .) will run 
through space values 1,2,3 only. 
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where are 5 linearly independent symmetric and traceless tensors, and E— is a multiple of the unit tensor S^j. 
S^ s \t) and <S( m )(t) are uniquely defined functions, whose explicit form depends on the particular representation of 
the _E-matrices chosen. A convenient one is the following: 



(S) 



E 










>': 


1 


i 










(2.16a) 




E 



(± 



1) = (l5_Y 
V32vr j 




E 



(±2) 
ij 




(2.16b) 



The excellence of this representation stems from its ability to display the spin features of the driving terms in (2.13) 
Such features are characterised by the relations 



, , * oo [p,<p) , E^ niUj = Y 2m {0, tp) 



(2.17) 



where n = x/|x| is the radial unit vector, and Yi m {6, (p) are spherical harmonics jL(|. Details about the a bove 
i?-matrices are given in Appendix A. In particular, the orthogonality relations (|A6|) can be used to invert (2.15): 



47T 

Y 

8tt 



4ps l3 (t) 
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-2, 



(2.18a) 
(2.18b) 



where an asterisk denotes complex conjugation. Note that S^ s \t) = \f^S{t) /3, where S{t) = SijSij(t) is the 
tensor's trace. 

We now take advantage of (2.15) to express the GW tidal force ( 2.1 3| ) as a sum of split terms like ( |2.2| ): 



f(x,<) = f( s )(x) 5 ( s )(;) 



(m) 



(2.19) 



with 



(S) 



g (s \t) = ^E^R 0m (t)ci 



(2.20a) 



(m) 
ij ' 



(2.20b) 



Straightforward application of (2.9) yields the formal solution of the antenna response to a GW perturbation: 

2 



N 



(2.21) 



with the n otatio n of ( |2.6| ) and (2.10) applied mutatis mutandi to the terms in ( 2.20| ). 

Equation ( 2.21 ) gives the response of an arbitrary elastic solid to an incoming weak GW, independently of the 
underlying gravity theory, be it General Relativity (GR) or indeed any other metric theory of the gravitational 
interaction. It is also valid for any antenna shape and any boundary conditions, thus giving the formalism, in particular, 
the capability of being used to study the response of a detector which is suspended by means of a mechanical device 
in the laboratory site — a situation of much practical importance. It is therefore very general. 

Equation ( p. 21 ) also tells us that that only monopole and quadrupole detector modes can possibly be excited by a 
metric GW. The nice thing about ( 2.21 ) is that it fully displays the monopole-quadrupole structure of the solution 
to our fundamental differential equations. 

In a non-symmetric body, all (or nearly all) the modes have monopole and quadrupole moments, and ( [2.21 ) precisely 
shows how much each of them contributes to the detector's response. A homogeneous spherical ante nna, which is 
very symmetric, has a set of vibrational eigenmodes which are particularly well matched to the form (2.21): it only 
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possesses one series of monopole modes and one (five-fold degenerate) series of quadrupole modes — see next section 
and Appendix B for details. The existence of so few modes which couple to GWs means that all the absorbed 
incoming radiation energy will be distributed amongst those few modes only, thereby making the sphere the most 
efficient detector, even from the sensitivity point of view. The higher energy cross section per unit mass repo rted 
for spheres on the basis of GR JlCfj , for example, finds here its qualitative explanation. The generality of (2.21), on 
the other hand, means that this excellence of the spherical detector is there independently of which is the correct GW 
theory. 

Before going further, let me mention another potentially useful application of the formalism so far. Cylindrical 
antennas, for in stance, are usually studied in the thin rod approximation; although this is generally quite satisfactory, 
equation (2.21) offers the possibility of eventually considering corrections to such simplifying hypothesis by use of 
more realistic eigenfunctions, such as those given in Recent new proposals for stumpy cylinder arrays p~9| ] 

may well benefit from the above approach, too. 



III. THE SPHERICAL ANTENNA 



To explore the consequences of (2.21) in a particular case, the mode amplitudes ujv(x) and frequencies lun must be 
specified. From now on I will focus on a homogeneous sphere whose surface is free of tractions and/or tensions; the 
latter happens to be quite a good approximation, even if the sphere is suspended in the static gravitational field poj . 

The normal modes of the free sphere fall into two families: so called toroidal — where the sphere only undergoes 
twistings which keep its shape unchanged throughout the volume — and spheroidal pjj ], where radial as well as 
tangential displacements take place. I use the notation 



nlm 



(x)e 



nlm 



(x) 



(3.1) 



for them, respectively; note that the index N of the previous section is a multiple index {nlm} for each family; I and 
m are the usual multipole indices, and n num bers from 1 to oo each of the /-pole modes. The frequencies happen to 
be independent of m, and so every one mode (3.1) is (2/+l)-fold degenerate. Further details about these eigenmodes 
are given in Appendix B. 



In order to see what (2.21) looks like in this case, integrals of the form (2.6) ought to be evaluated. It is straightfor- 
ward to prove that they all vanish for the toroidal modes, the spheroidal modes contributing the only non-vanishing 
terms; after some algebra one finds 



where 



AS) 

J nlm 



Am') 
•I nlm 



l 

M 
1 

M 



u^^(x) • f (s) (x) d 3 x = a n 5io S m0 



Sphere 



f u r ^„(x).f( m ')(x)d 3 x = &„ 

J Sphere 



$12 $n 



(3.2a) 



(3.2b) 



(3.3a) 



1 r R 

-— J [A n2 (r) + 3 B n2 (r)} pr 3 dr 



(3.3b) 



The functions A„i(r), B„i(r) are given in Appendix B, and R is the sphere's radius. To our reas suran ce, o nly th e 
monopole and quadrupole sphere modes survive, as seen by the presence of the factors Siq and 812 in ( 3.2a ) and ( 3.21: ), 
respectively. The final series is thus a relatively simple one, even in spite of its generalityn: 



u 0M) = y] — u noo(x).gio ) W + Y] 

— tOnO 



Un2 



Y U «2m(x) gffl (t) 



,m=-2 



(f >0) 



(3.4) 



From now on I will drop the label P, meaning spheroidal mode, to ease the notation since toroidal modes no longer appear 
in the formulae. 
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where, it is recalled, 



lf' m) (i)= / g (s ' m \t') sincj nl (t-t')dt' , (m=-2,...,2) (3.5) 



9nl 

Equation ( |3.4| ) constitutes the sphere's response to an arbitrary tidal GW perturbation, and will be used to analyse 
the sensitivity of the spherical detector in the next section. Before doing so, however, a few comments on the antenna's 
signal deconvolution capabilities, within the context of a completely general metric theory of GWs, are in order. 



A. The deconvolution problem 



Let us first of all take the Fourier transform of (p\4|): 

/oo 
u(x, t) e~ lut dt 
-OO 

This is seen to be 



(3.6) 



OO 

U(x,cj) = - V" — u n0 o(x)G (s) (w) [8(cj-u n o) - S(u)+u) n0 )] 

7 * ■* 



I ^ UJ n0 
n—1 



7T \ - b n 



n=l 



0Jn2 



2 

E 

.m=-2 



u„ 2m (x)G^M 



[S(oj-u> n2 ) - 5{uj+uj n2 )] 



where G^(uj) and G^ m \uj) are the Fourier transforms of g^ s \t) and g( m \t), respectively: 



G (s < m) (w) 



'(*)< 



dt 



(3.7) 



(3.8) 



The 5-function factors are of course idealisations corresponding to infinitely long integration times and infinitely 
narrow resonance linewidths — but the essentials of the ensuing discussion will not be affected by those idealisations. 

If the measuring system were (ideally) sensitive to all frequencies, filters could be applied to examine the antenna's 
oscillations at each monopole and quadrupole frequency: a single transducer would suffice to reveal G^ s ' (w) around the 
monopole frequencies ui n0 , whilst five (placed at suitable positions) would be required to calculate the five degenerate 
amplitudes G^ m \io) around the quadrupole frequencies w„2- Once the six functions G^ s,m \uj) would have thus been 
determined, inverse Fourier transforms would give us the func tions g^ s (t) , and thereby the six Riemann tensor 
components Roioj(t) through inversion of the second equations ( 2.20| ), i.e., as an expansion like (2.15) — only with g's 
instead of 6>'s. Deconvolution would then be complete. 

Well, not quite. . . Knowledge of the Riemann tensor in the laboratory frame coordinates is not really sufficient to 
say the waveform has been completely deconvolved, unless we also know the source position in the sky. There clearly 
are two possibilities: 

i) The source position is known ahead of time by some other astronomical observation methods. Let me rush to 
emphasise that, far from trivial or uninteresting, this is a very important case to consider, specially during the 
first stages of GW Astronomy, when any reported GW event will have to be thoroughly checked by all possible 
means. 

If the incidence direction is known, then a rotation must be applied to the just obtained quantities Roioj(t), 
which takes the laboratory z-axis into coincidence with the incoming wave propagation vector. A classification 
procedure must thereafter be applied to the so transformed Riemann tensor in order to see which is the theory 
(or class of theories) compatible with the actual observations. Such classification procedure has been described 
in detail in [^2| . 

The spherical antenna is thus seen to have the capability of furnishing the analyst sufficient information to 
discern amongst different competing theories of GW physics, whenever the wave incidence direction is known 
prior to detection. 



G 



ii) The source position is not known at detection time. This makes things more complex, since the above rotation 
between the laboratory and GW frames cannot be performed. 

In order to deconvolve the incidence direction in this case, a specific theory of the GWs must be assumed — a 
given choice being made on the basis of whatever prior information is available or, simply, dictated by the the 
decision to probe a particular theory. Wagoner and Paik |?J propose a method which is useful both for GR 
and BD theory, their idea being simple and elegant at the same time: since neither of these theories predicts 
the excitation of the m=±l quadrupole modes of the wave, the source position is determined precisely by the 
rotation angles which, when applied to the laboratory axes, cause the amplitudes of those antenna modes to 
vanish; the rotated frame is thereby associated to the GW natural frame. 

A generalisation of this idea can conceivably be found on the basis of a detailed — and possibly rather casuistic — 
analysis of the canonical forms of of the Riemann tensor for a list of theories of gravity, along the following line 
of argument: any one particular theory will be characterised by certain (homogeneous) canonical relationships 
amongst the monopole and quadrupole components of the Riemann tensor, g^ S ' m ^(t), and so enforcement of 
those relations upon rotation of the laboratory frame axes should enable determination of the rotation angles 
or, equivalently, of the incoming radiation incidence direction. Scalar-tensor theories e.g. have g^ 1 '^) = in 
their canonical forms, hence Wagoner and Paik's proposal for this particular case. 

Before any deconvolution procedure is triggered off, however, it is very important to make sure that it will be 
viable. More precisely, since the transformation from the laboratory to the ultimate canonical frame is going 
to be linear, invariants must be preserved. This means that, even if the source position is unknown, certain 
theories will forthrightly be vetoed by the observed Roioj (t) if their predicted invariants are incompatible with 
the observed ones. To give but an easy example, if RoiOj(t) is observed to have a non-null trace Rom(t), then a 
veto on GR will be readily served, and therefore no algorithm based on that theory should be applied. 

I would like to make a final remark here. Assume a direction deconvolution procedure has been successfully 
carried through to the end on the basis of certain GW theory, so that the analyst comes up with a pair of 
numbers (9, tp) expressing the source's coordinates in the sky. Of course, these numbers will represent the 
actual source position only if the assumed theory is correct. Now, how do we know it is correct? Strictly 
speaking, "correctness" of a scientific theory is an asymptotic concept — in the sense that the possibility always 
remains open that new facts be eventually discovered which contradict the theory — , and so reliability of the 
estimate (9, tp) of the source position can only be assessed in practice in terms of the consistency between the 
assumed theory and whatever experimental evidence is available to date, including, indeed, GW measurements 
themselves. It is thus very important to have a method to verify that the estimate {9, ip) does not contradict 
the theory which enabled its very determination. 

Such verification is a logical absurdity if only one measurement of position is available; this happens for instance 
if the recorded signal is a short burst of radiation, and so two antennas are at least necessary to check consistency 
in that case. The test would proceed as a check that the time delay between reception of the signal at both 
detectors is consistent with the calculated (9, (p)^ given their relative position and the wave propagation speed 
predicted by the assumed theory. If, on the other hand, the signal being tracked is a long duration signal, then 
a single antenna may be sufficient to peform the test by looking at the observed Doppler patterns and checking 
them against those expected with the given (9, tp). 

The above considerations have been made ignoring noise in the detector and monitor systems. A fundamental 
constraint introduced by noise is that it makes the antenna bandwidth limited in sensitivity. As a consequence, any 
deconvolution procedure is deemed to be incomplete or, rather, ambiguous p3t], since information about the signal 
can possibly be retrieved only within a reduced bandwidth, whilst the rest will be lost. I thus come to a detailed 
discussion of the sensitivity of the spherical GW antenna in the next section. 

IV. THE SENSITIVITY PARAMETERS 

I will consider successively amplitude and energy sensitivities; the first leads to the concept of transfer function, 
while the second to that of absorption cross section. I devote separate subsections to analyse each of them in some 
detail. 



3 Note that the two detectors will agree on the same (9, tp), even if the assumed theory is wrong, since the sphere deformations 
will be the same if caused by the same signal. 
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A. The transfer function 



A widely used and useful concept in linear system theory is that of transfer function p3| . It is defined as the Fourier 
transform of the system's impulse response, or as the system's impedance/admittance, and can be inferred from the 
frequency response function (3.7). 



We recall from the previous section that the sphere is a multimode device — due to its monopole and five-fold 
degenerate quadrupole modes. It appears expedient to define a multimode or vector transfer function as a useful 
construct which encompasses all six different modes into a single conceptual block, according to 



U(x, uj)=J2 z(Q) ( x > ^) G(a) 



(4.1) 



where G^(to) are the six driving terms G^ S ' m \uj) given in ( |3.8|) . The transfer function is Z'") (x, w), and its 
"vector" character alluded above is reflected by the multimode index a. Looking at (3.7) it is readily seen that 



Z (5) (x,u) = — } — u„oo(x) [S(cu-LOno) - S(uj+uj n0 )} 



71=1 



(4.2a) 



7 

Z< m )(x, w ) = - V -=-u„ 2m (x) Mw-wna) - 6(u,+u; n2 )} (m = -2,...,2) 



n=l 



(4.2b) 



As we observe in these formulae, the sphere's sensitivity to monopole excitations is governed by a n /u> n o, and to 
quadrupole ones by b n /u> n 2- Closed expressions happen to exist for a n and b n ; using the notation of Appendix B, 
they are 



On _ 3C(n, 0) j 2 (gn0^) 
R 8tt q n0 R 



(4.3a) 



bn 

R 



3C*(n,2) 
8^ 



/3 3 (fe„ 2 -R) 5 3 -j — f3 1 {q n2 R) — — — 

9n2-K K„2 K n 2K 



(4.3b) 



Numerical investigation of the behaviour of these coefficients shows that they decay asymptotically as n 2 : 



, b„ n — const x n 2 



(4.4) 



Likewise, it is found that the frequencies io n o and uj n 2 diverge like n for large n, so that Z'°'(x,w) drops as uj~ 3 for 
large w. Figures 6 and 7 display a symbolic plot of uj 3 (x, uj) and w 3 Z' m ^ (x, lj) , respectively, which illustrates the 
situation: monopole modes soon reach the asymptotic regime, while there appear to be 3 subfamilies of quadrupole 
modes regularly intertwined; the asymptotic regime for these subfamilies is more irregularly reached. Note also the 
perfectly regular alternate changes of phase (by tt radians) in both monopole and each quadrupole family. 

The sharp fall in sensitivity of a sphere for higher frequency modes (n~ 3 ) indicates that only the lowest ones stand 
a chance of being obervable in an actual GW antenna. I report in Table I the numerical values of the relevant 
parameters for the first few monopole and quadrupole modes. The reason for the last (fourth) columns will become 
clear later. 



B. The absorption cross section 



Let us calculate now the energy of the oscillating sphere. We first define the spectral energy density at frequency 
u>, which is naturally given byn 



W(uj) 



1 
T 



\j |U(x,^)| 2 pd 3 x 

Solid * 



(4.5) 



T is the integration time — assumed very large. The peaks in the ^-functions diverge like T/tt. 
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TABLE I. First few monopole (left) and quadrupole (right) sphere parameters, for a <r=0.33 material. First and second 
columns on either side of the central line number the modes and give the corresponding eigenvalue; rows are intertwined in 
order of ascending fre quen cy, whi ch is proportional to kR — see ( JBq ) below. Third columns contain the a n and b n coefficients 
defined in equations (3.3a) and (3.3b), respectively; the fourth columns display the cross section ratios (kioai/k„oa n ) 2 and 



(&12&1 / k„2b n ) for higher frequency modes, respectively, taking as reference the lowest in each familiy — cf. equations (4.f8) 



n 


k„oR 


a„/R 


Clo/CnO 


n 


k„2R 


b n /R 


0"l2/cr n 2 










1 


2.650 


0.328 


1 










2 


5.088 


0.106 


2.61 


1 


5.432 


0.214 


1 


















3 


8.617 


-1.907xl0~ 2 


27.95 










4 


10.917 


-9.101xl0" 3 


76.42 


2 


12.138 


-3.772xl(T 2 


6.46 


















5 


12.280 


1.387xl0~ 2 


25.99 










6 


15.347 


6.879xl0~ 3 


67.87 


3 


18.492 


1.600 xl(T 2 


15.49 











and can be easily evaluated: 

oo 2 
W(cj) = -ttM < GW(u) [6(u;-u n0 ) + S(u+oj n0 )] + 



1 OO 



|G (m) (w) 



[d{uj-uj n2 ) + 6{u)+u) n2 )} 



(4.6) 



The energy at any one spectral frequency tu n i is obtained by integration of the spectral density in a narrow interval 
around u) = ±uJ n r- 



E(u> n l) 



+ / W{u>) ^ 



2TT 



In particular, 



E(Lo n0 ) = ~Mc& |G( S) Ko) 



E{u n2 ) = \Mb 2 n £ |g^K 2 ) 



m=-2 



(4.7) 

(4.8a) 
(4.8b) 



The sensitivity parameter associated with the vibrational energy of the modes is the detector's absorption cross 
section, denned as the energy it absorbs per unit incident GW spectral flux density, or 



$(cj) 



(4.9) 



where $(w) is the number of joules per square metre and Hz carried by the GW at frequency lo as it passes by the 
antenna. Thus, for the frequencies of interest, 



o"ab s (^«o) = \ Ma n 



|G (S) Ko) 



i lfL2 EL-2 |g(" 1 W)| 



(4.10a) 



(4.10b) 



These quantities have very precise values, but such values can only be calculated on the basis of a specific underlying 
theory of the GW physics. In the absence of such theory, neither <&(cj) nor G^ s - m \uj) can possibly be calculated, since 
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they are not theory independent quantities. To date, only GR c alcula tions have been reported in the literature ]7|,p|,^0| 



As I will now show, even though the fractions in the rhs of ( 4.10 ) are not theory independent, some very general 
results can still be obtained about the sphere's cross section within the context of metric theories of the gravitational 
interaction. To do so, it will be necessary to go into a short digression on the general nature of weak metric GWs. 

No matter which is the (metric) theory which happens to be the "correct one" to describe gravitation, it is beyond 
reasonable doubt that any GWs reaching the Earth ought to be very weak. The linear approximation should therefore 
be an extremely good one to describe the propagating field variables in the neighbourhood of the detector. In such 
circumstances, the field equations can be derived from a Poincare invariant variational principle based on an action 
integral of the type 

C{ijj A ,ip Alll )d 4 x (4.11) 

where the Lagrangian density £ is a quadratic functional of the field variables ^a (x) and their space-time derivatives 
ipA,n( x )'i these variables include the metric perturbations h^, plus any other fields required by the specific theory 
under consideration — e.g. a scalar field in the theory of Brans-Dicke, etc. The requirement that C be quadratic 
ensures that the Euler-Lagrange equations of motion are linear. 

The energy and momentum transported by the waves can be calculated in this formalism in terms of the components 
t^ v of the canonical energy-momentum tensor^ 

BC 

^ V ^)=Hm — ^a-^ V (4-12) 

The flux energy density, or Poynting, vector is given by Sj — c 2 r° l , i.e., 

BC 

S(x,i)= C 3 ]T -pVVu (4.13) 



A 



A 



where "= d/dt. Any GW hitting the antenna will be seen plane, due to the enormous distance to the source. If k 
is the incidence direction (normal to the wave front), then the fields will depend on the variable ct— k-x, so that the 
GW energy reaching the detector per unit time and area is 

BC 

0(t)^k-S(x,t) = -c 2 J2 -77- (4.14) 

A d ^ A 

where x is the sphere's centre position relative to the source — which is fixed, and so its dependence can be safely 



dropped in the lhs of the above expression. The important thing to note in equation (4.14) is that it tells us that 4>(t) 
can be written as a quadratic form in the time derivatives of the fields ipA- As a consequence, the spectral density 
&(u>), defined by 



f°° _ , ,dw 



4>{t)dt= / $M— (4.15) 

J — oo 

can be ascertained to factorise as 

= lu 2 $ (w) (4.16) 

where $q{lo) is again a q uadr atic function of the Fourier transforms ^a(^) of the fields 4>A- O n the other hand, 
the functions G^ S ' m \uj) in ( 4.10 ) which, it is recalled, are the Fourier transforms of g( s,m \t) in ( [2.20 ), contain second 



order derivatives of the metric fields hp U , and therefore of all the fields ipA as a result of the theory's field equations. 
Since we are considering plane wave solutions to those equations, all derivatives can be reduced to time derivatives 
— just like in ( 4.14| ) above. We can thus write 



5 This tensor is not symmetric in general, but can be symmetrized by a standard method due to Belinfante |[25|]26| . For 
the considerations which follow in this paper it is unnecessary to go into those details, and the canonical form ( ETTgM will be 
sufficient. 
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(4.17) 



with \E'( s, ' m )((jj) suitable linear combinations of the ^^(w). Replacing the last two equations into ( 4.1 (j ) and 
manipulating dimensions expediently, we come to the remarkable result that 



GMv 2 , 2 
^- L (k n oa n ) 



GMv 2 

CTabs(Wn2) = Kq(#) (fc n2 &n) 2 



(4.18a) 
(4.18b) 



where v\ = (2+2(7 ) 



being the speed of sound in the detector's material, and a its Poisson ratio; G is the 



Gravitational constant. The "remarkable" about the above is that the coefficients ifs(H) and Kq(H) are independent 
of frequency: they exclusively depend on the underlying gravitation theory, which I symbolically denote by H. To see 
that this is the case, it is enough to consider a monochromatic incident wave: since the coefficients Ks(&) and Kq(N) 
happen to be invariant with respect to field amplitude scalings, this means they will only depend on the amplitudes' 
relative weights, i.e., on the field equations' specific structure. 

By way of example, it is interesting to see what the results for General Relativity (GR) and Brans-Dicke (BD) 
theory are. After somehow lengthy algebra it is found that 



N = GR 



and 



K S (X) = 



16 7T 2 

15 



(4.19) 



^s(K) = ^-(3 + 2ft)- 



K = BD 



kit 



(3+2C) 2 



(4.20) 



jr g (N) = !ff- [l + i(3 + 2f7)- 2 fc] 



kn 



{3+2n) 2 



In the latter formulae, fl is the usual Brans-Dicke parameter oj f27fl , renamed here to avoid confusion with frequency, 
and k is a dimensionless parameter, generallyof order one, depending on the source's properties pSfl . As is well known, 
GR is obtained in the limit Q — > oo of BD fll5|] ; the quoted results are of course in agreement with that limit. 

Incidentally, an interesting consequence of the above equations is this: though not explicitly shown in this paper 
(see, however, reference 0), the presence of a scalar field in the theory of Brans and Dicke causes not only the 
monopole sphere's modes to be excited, but also the m=0 quadrupole ones; what we see in equations ( 4.20[ ) is that 
precisely 5/6 of the total energy extracted from the scalar wave goes into the antenna's monopole modes, whilst there 
is still a remaining 1/6 which is communicated to the quadrupoles, independently of the values of and Aj^j. This 
somehow non-intuitive result finds its explanation in the structure of the Riemann tensor in BD theory, in which 

(S) 

the excess Rowj with respect to General Relativity happens not to be proportional to the scalar part E- , but to a 
combination of EST 1 and E$ . 



Equations ( 4.18 ) show that, no matter which is the gravity theory assumed, the sphere's absorption cross sections 
for higher modes scale as the successive coefficients (k n oa n ) 2 and (k n 2b n ) 2 for monopole and quadrupole modes, 
respectively. In particular, the result quoted in [fiof that cross section for the second quadrupole mode is 2.61 times 
less than that for the first, assuming GR, is in fact valid, as we now see, independently of which is the (metric) theory 
of gravity actually governing GW physics. The fourth columns in Table I display these scaling properties. It is seen 
that the drop in cross section from the first to the second monopole mode is as high as 6.46. It should however be 
stressed that the frequency of such mode would be over 4 kHz for a (likely) sphere whose fundamental quadrupole 
frequency be 900 Hz [Q. Note finally the asymptotic cross section drop as n~ 2 for large n — cf. equation (4.4) and 
the ensuing paragraph. 



Note however that since monopole and quadrupole detector modes occur at different frequencies, this particular distribution 
of energy may not be seen if the sphere's vibrations are monitored at a single resonance. 
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V. CONCLUSION 



The main purpose of this paper has been to set up a sound mathematical formalism to tackle with as much 
generality as possible any questions related to the interaction between a resonant antenna and a weak incoming GW, 
with much special emphasis on the homogeneous sphere. New results have been found along this line, such as the 
scaling properties of cross sections for higher frequency modes, or the sensitivity of the antenna to arbitrary metric 
GWs; also, new ideas have been put forward regarding the direction deconvolution problem within the context of 
an arbitrary metric theory of GW physics. Less spectacularly, the full machinery has also been applied to produce 
independent checks of previously published results. 

The whole investigation reported herein has been developed with no a priori assumptions about any specific 
(metric) theory of the GWs, and is therefore very general. "Too general solutions" are often impractical in science; 
here, however, the "very general" appears to be rather "cheap", as seen in the results expressed by the equations 
of section 3 above. An immediate consequence is that solid elastic detectors of GWs (and, in particular, spheres) 
offer, as a matter of principle, the possibility of probing any given theory of GW physics with just as much effort 
as it would take, e.g., to probe General Relativity: the vector transfer function of section 4 supplies the requisite 
theoretical vehicle for the purpose. 

An important question, however, has not been considered in this paper. This is the transducer problem: the sphere's 
oscillations can only be revelaed to the observer by means of suitable (usually electromechanical) transducers. These 
devices, however, are not neutral, i.e., they couple to the antenna's motions, thereby excercising a back action on it 
which must be taken into consideration if one is to correctly interpret the system's readout. Preliminary studies and 
proposals have already been published but further work is clearly needed for a more thorough understanding of 
the problems involved. 

Progress in this direction is currently being made — which I expect to report on shortly. The formalism developed 
in this paper provides basic support to that further work. 
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APPENDIX A: 



Let e x ,e y ,e z be three orthonormal Cartesian vectors defining the sphere's laboratory reference frame. We define 
the equivalent triad 



e<°> = e z , - -i 



-^(e x ±ie y ) (Al) 



having the properties 



e *(m') . e (m) = ( m>m ' = -1,0,1 (A2) 

We say that the vectors (Al) are the natural basis for the 1=1 irreducible representation of the rotation group; they 



behave under arbitrary rotations precisely like the spherical harmonics Yi m (6, ip). In particular, if a rotation of angle 
a around the z-axis is applied to the original frame, then 



e 



( ±1 )^exp(±m)e( ±1 ) , e<°>-e<°> (A3) 



Higher rank tensors have specific multipole characteristics depending on the number of tensor indices, and the 
above basis lends itself to reveal those characteristics, too. For example, the five dimensional linear space of traceless 
symmetric tensors supports the 1=2 irreducible representation of the rotation group, while a tensor's trace is an 
invariant. A general symmetric tensor can be expressed as an "orthogonal" sum of a traceless symmetric tensor and 
a multiple of the unit tensor. A convenient basis to expand any such tensor is the following: 
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,(0) 



,(<>) 



e (0) ^e'- 1 ' +e(- 1 ) 



,(o) 



2 (i) (8e (- 1 )+e(- 1 )®e( 1 )-2e(°) 



,(i) 



3 (D 



+ e 



(0) 



ie (0) 
P (0) 



(A4a) 
(A4b) 
(A4c) 
(A4d) 



The elements (A4a) get multiplied by e ±2lc * in a rotation of angle a around the 2-axis, respectively, the ( A4b| ) by 



e , and (A4c) and (A4d) are invariant, as is read ily seen . These properties define the "spi n characteristics" of the 



corresponding tensors. Also, the five elements (A4a)-(A4c) are traceles s te nsors, while (A4d) is the unit tensor. Any 
symmetric tensor can be expressed as a linear combination of the six (A4), and the respective coefficients carry the 
information a bout the weights of the different monopole and quadrupole components of the tensor. 

Equations (2. IE) in the text are the matrix representation of the above tensors in the Cartesian basis e x ,e y ,e z , 
except that they are multiplied by suitable coefficients to ensure that the conditions 



E iP n i n 3 = Y oo(0, <p) , E^' mrij = Y 2m (9, (p) 



(A5) 



where n = x/|x| is the radial unit vector, hold. They are arbitrary, but expedient for the calculations in this paper. 
The following orthogonality relations can be easily established: 



15 r p (S) p (m) 

8tt ° m m ' ^ij C 'J 



E (S) E iS) 

1 V V 



4tt 



(A6) 



with the indices m,m' running from —2 to 2, and with an understood sum over the repeated i and j. It is also easy 
to prove the closure properties 



E %] E kl 



m=— 2 



*(m) n,(m) 



E kl' = 8^ i^Sjl + SiiSjk) 



(A7) 



Equations (A6) and (A7) constitute the completeness equations of the i^-matrix basis of Euclidean symmetric 
tensors. 



APPENDIX B: 



This Appendix is intended to give a rather complete summary of the frequency spectrum and eigenmodes of a 
uniform elastic sphere. Although this is a classical problem in Elasticity Theory ^9|, some of the results which follow 
have never been published so far. Also, its scope is to serve as reference for notation, etc., in future work. 

The uniform^ elastic sphere's normal modes are obtained as the solutions to the eigenvalue equation 



/xV 2 u + (A + n) V(V-u) = -uj 2 qu 



(Bl) 



with the boundary conditions that its surface be free of any tensions and/or tractions; this is expressed by the 
equations Q 



Uij nj = at r = R 
where R is the sphere's radius, n the outward normal, and cr^ the stress tensor 

a l j = Xu kk S i: j + 2/j.Uij 



(B2) 



(B3) 



with = ^{uij + Uj i), the strain tensor, and A, [i the Lame coefficients pJ| . 

Like any differentiable vector field, u(x) can be expressed as a sum of an irrotational vector and a divergence-free 
vector, 



u(x) = u irrot .(x) + u div _f ree (x) , 



(B4) 



7 By uniform I mean its density g is constant throughout the solid in the unperturbed state. 
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say; on substituting this into equation (Bl), and after a few easy manipulations, one can sec that 

(V 2 + k 2 ) u div _ frec (x) = , (V 2 + q 2 ) u irrot . (x) = 

where 



QLO 



q 



QUO 



\ + 2fi 

Now the irrotational component can generically be expressed as the gradient of a scalar function, i.e., 

Uirrot. (x) = V0(x) 

while there are two linearly independent divergence-free components which, as can be readily verified, are 



u 



(1) 



- free 



(x)=L^W(x), and (x) = V x LV>< 2 >(x 



(B5) 

(B6) 

(B7) 
(B8) 



wher e L = — ixx V is the "angular momentum" operator, cf. |T(|, and ip^ and ip^ are also scalar functions. If (B7) 
and (B8) are now respectively substituted in (B5), it is found that </>(x), ^^(x), and ^^(x) satisfy Hclmholtz 
equations: 



(V 2 + fc 2 ) ?/>(x) = , (V 2 + q 2 ) 0(x) = 
where stands for either ^W(x) or ?/>( 2 )(x). Therefore 

</>(x) = ji{qr)Y [m (n) , ^>(x) = j t (kr) Y lm (n) 



(B9) 



(BIO) 



in order to ensure regularity at the centre of the sphere, r=0. Here, j\ is a spherical Bessel function — see p0| for 
general conventions on these functions — , and Yi m a spherical harmonic fig]. Finally thus, 



u(x) = ^ V0(x) + ^ L^(x) +^Vx L, -(xi 



(Bll) 



where Co,Ci,C 2 are three constants which will be determined by the boundary conditions (|B2| ) (the denominators 
under them have been included for notational convenience). After lengthy algebra, those conditions can be expressed 
as the following system of linear equations: 



fo{qR) ~ ^ q 2 R 2 P (qR)] C - 1(1 + 1) Pi{kR) C 2 = 



0i(qR)C o 



where 



lf3 2 {kR) + {^--l} fo(kR) 



(3 1 (kR)C 1 = 
C 2 = 



d 

dz 



dz 2 



(B12a) 
(B12b) 
(B12c) 

(B13) 



There are clearly two families of solutions to ( B12 ): 
i) Toroidal modes. These are characterised by 

0i(kR) = O, C =C 2 =0 



(B14) 



The frequencies of these modes are independent of A, and thence independent of the material's Poisson ratio. 
Their amplitudes are 



with 



u n; m ( x ) = T m(r) iLY lm (n) 



T n i(r) = Cx(n,l)ji(k n ir) 



(B15) 



(B16) 



and Ci(n,l) a dimensionl ess n ormalisation constant determined by the general formula ( J2.8| ); k n iR is the n-th 
root of the first equation (B14) for a given I. 
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ii) Spheroidal modes. These correspond to 



dct 



Pi(qR) 



1(1 + 1) 0x(kR) 
\f3 2 (kR) + {^ 1 -l) o (kR) 



= 



and C\ = 0. The frequencies of these modes do depend on the Poisson ratio, and their amplitudes are 

u ni m ( x ) = A m (r) Yim (n) n - B n i (r) in x LY lm (n) 
where A n i(r) and B n i(r) have the somewhat complicated form 

A nl (r) = C(n,l) 



mmR)3'Mnir) - 1(1 + 1) HlmR) JJ ^- 



B„i(r) = C(n,l) 



a jj\ Ji(qmr) q n i R ( m {Kir ji(k n ir)} 

P3(kniR) — Pi(q n iR) : 

q n ir k n i k n ir 



with accents denoting derivatives with respect to implied (dimensionless) arguments, 



(BIT) 

(B18) 

(B19a) 
(B19b) 

(B20) 



and C(n,l) a ne w nor malisation constant. It is understood that q n i and k n i are obtained after the (transcen- 
dental) equation (B17) has been solved for u> — cf. equation (B6). 



In actual practice equations (B14) and (B17) are solved for the dimensionless quantity kR, which will hereafter be 
called the eigenvalue. In view of (B6), the relationship between the latter and the measurable frequencies (in Hz) is 
given by 



2tt 



Q R 2 



1/2 



kR 
2^ 



(B21) 



It is more useful to express the frequencies in terms of the Poisson ratio, er, and of the speed of sound v s in the 
selected material. For this the following formulas are required — see e.g. pJ| : 

(B22) 



where Y is the Young modulus, related to the Lame coefficients and the Poisson ratio by 



Y 



(3A + 2fj) n 



Hence, 



2(l + cr)^ 



(kR) 



X 



2(A + M ) 



2ttvT 



R 



(B23) 



(B24) 



Equation ( B24 ) provides a suitable transformation formula from abstract number eigenvalues (kR) into physical 
frequencies is, for given material's properties and sizes. 

Tables [H and [n| respectively display a set of values of (kR) for toroidal and spheroidal modes. While GWs can only 
couple to quadrupolc and monopole modes, it is important to have some detailed knowledge of analytical results, as the 
sphere's frequency spectrum is rather involved. It often happens, both in numerical simulations and in experimental 
determinations, that it is very difficult to disentangle the wealth of observed frequency lines, and to correctly associate 
them with the corresponding eigenmode. Complications are further enhanced by partial degeneracy lifting found in 
practice (due to broken symmetries), which result in even more frequency lines in the spectrum. Accurate analytic 
results should therefore be very helpful to assist in frequency identification tasks. 
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TABLE II. List of a few spheroidal eigenvalues, ordered in columns of ascending harmonics for each multipole value. 
Spheroidal eigenvalues depend on the sphere's material Poisson ratio — although this dependence is weak. In this table, 
values are given for a = 0.33. Note that the table contains all eigenvalues less than or equal to 11.024 yet is not exhaustive for 
values larger than that one; this would require to stretch the table horizontally beyond I — 10 — see Figure I for a qualitative 
inspection of trends in eigenvalue progressions. 



n 


I = 


/ = 1 


/ = 2 


I = 3 


I = 4 


I = 5 


I = 6 


I = 7 


I = 8 


I = 9 


I = 10 


1 


5.4322 


3.5895 


2.6497 


3.9489 


5.0662 


6.1118 


7.1223 


8.1129 


9.0909 


10.061 


11.024 


2 


12.138 


7.2306 


5.0878 


6.6959 


8.2994 


9.8529 


11.340 


12.757 


14.111 


15.410 


16.665 


3 


18.492 


8.4906 


8.6168 


9.9720 


11.324 


12.686 


14.066 


15.462 


16.867 


18.272 


19.664 


4 


24.785 


10.728 


10.917 


12.900 


14.467 


15.879 


17.243 


18.589 


19.930 


21.272 


22.619 


5 


31.055 


13.882 


12.280 


14.073 


16.125 


18.159 


19.997 


21.594 


23.043 


24.426 


25.778 



TABLE III. List of a few toroidal eigenvalues, ordered in columns of ascending harmonics for each multipole value. Unlike 
spheroidal eigenvalues, toroidal eigenvalues are independent of the sphere's material Poisson ratio. Note that the table contains 
all eigenvalues less than or equal to 12.866 yet is not exhaustive for values larger than that one; this would require to stretch 
the table horizontally beyond I = 11 — see Figure II for a qualitative inspection of trends in eigenvalue progressions. 



n 


I = 1 


/ = 2 


/ = 3 


I = 4 


I = 5 


I = 6 


I = 7 


I = 8 


I = 9 


I = 10 


I = 11 


1 


5.7635 


2.5011 


3.8647 


5.0946 


6.2658 


7.4026 


8.599 


9.6210 


10.711 


11.792 


12.866 


2 


9.0950 


7.1360 


8.4449 


9.7125 


10.951 


12.166 


13.365 


14.548 


15.720 


16.882 


18.035 


3 


12.323 


10.515 


11.882 


13.211 


14.511 


15.788 


17.045 


18.287 


19.515 


20.731 


21.937 


4 


15.515 


13.772 


15.175 


16.544 


17.886 


19.204 


20.503 


21.786 


23.055 


24.310 


25.555 


5 


18.689 


16.983 


18.412 


19.809 


21.181 


22.530 


23.860 


25.174 


26.473 


27.760 


29.035 



In Figures 1 and 2 a symbolic line diagramme of the two families of frequencies of the sphere's spectrum is presented. 
Spheroidal eigenvalues have been plotted for the Poisson ratio u= 0.33. Although only the 1=0 and 1=2 spheroidal 
series couple to GW tidal forces, the plots include other eigenvalues, as they can be useful both in bench experiments 
— cf. equation ( 2.12| ) above — and for vetoing purposes in a spherical antenna. 

Figures 3, 4 and 5 contain plots of the first three monopole and quadrupole functions T„;(r), A n i(r) and B n i(r), 
always for <t=0.33. T n o(r) and B„o(r) have however been omitted; this is because they are multiplied by an identically 
zero angular coefficient in the amplitude formulae ( B15| ) and (B18). Indeed, monopole vibrations are spherically 
symmetric, i.e., purely radial. 
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List of Figures 



Figure 1 The homogeneous sphere spheroidal eigenvalues for a few multipole families. Only the 1=0 and 1=2 
families couple to metric GWs, so the rest are given for completeness and non-directly-GW uses. Note that there 
are fewer monopole than any other Z-pole modes. The lowest frequency is the first quadrupole. The diagramme 
corresponds to a sphere with Poisson ratio <7=0.33. Frequencies can be obtained from the plotted values through 
equation (B6) for any specific case. 



Figure 2 The homogeneous sphere toroidal eigenvalues. None of these couple to GWs, but knowledge of them can 
be useful for vetoing purposes. Thes e ei genvalues are independent of the material's Poisson ratio. To obtain actual 
frequencies from plotted values, use (B6). The lowest toroidal eigenvalue is kR = 2.5011, with 1=2, and happens to 
be the absolute minimum sphere's eigenvalue. Compared to the spheroidal kR = 2.6497, also with 1=2, its frequency 
is 5.61% smaller. Note also that there are no monopole toroidal modes. 



Figure 3 First three spheroidal monopole radial functions A n o(r) (n = 1,2,3), equation ( B19a ) 



Figure 4 First thr ee sp heroidal quadrupole radial functions A n 2(r) (continuous line) and B n 2(r) (broken line) 
(n = 1,2,3), equations (B19). 



Figure 5 First three toroidal quadrupole radial functions T„ 2 (r) (n = 1, 2, 3), equation ( B16| ). A common feature 
to these radial functions (also in the two previous Figures) is that they present a nodal point at the origin (7' = 0), 
while the sphere's surface (r/R = 1) has a non-zero amplitude value, which is largest (in absolute value) for the lowest 
n in each group. 



Figure 6 The scalar component 7S S ^ (x, uj) of the multimode transfer function, ( 4.2a ). The diagramme actually 
displays uj 3 Z( s )(x, uj), so asymptotic behaviours are better appreciated. It is given in units of /_t/ pR, and a factor 
(ir/i) u n oo(x), the eigenmode amplitude, has been o mitt ed, too. <5-function amplitudes are symbolically taken as 1. 
Note that the asymptotic regime, given by equation (4.4), is quickly reached. 



Figure 7 The quadrupole component Z^ m ^(x, uj) of the multimode transfer function, (4.2b). The same prescriptions 
of Figure 6 apply here; the plot is therefore independent of the value of m. Note the presence of three subfamilies of 
peaks; asymptotic regimes are reached with variable speed for these subfamilies, and less rapidly than for monopole 
modes, anyway. 
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FIG. VI 
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Quadrupole modes 
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